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We show that the electronic structures of the title compounds predicted by density functional 
theory (DFT) are well described by tight binding models. We determine the frustration ratio, J' / J, 
of the Heisenberg model on the anisotropic triangular lattice, which describes the spin degrees of 
freedom in the Mott insulating phase for a range of Pd(dmit)2 salts. All of the antiferromagnetic 
materials studied have J' / J < 0.5 or J' / J > 0.9 and all salts with 0.5 < J' / J < 0.9 are known, 
experimentally, to be charge ordered, valence bond solids or spin liquids. 



The interplay of geometrical frustration and elec- 
tronic correlations produces a wide range of exotic 
phenomena [1, 2] in the organic charge transfer salts 
Me4-nEtnPn[Pd(dmit)2]2 (henceforth Pn-n) [3]. At 
ambient pressure and low temperature these materials 
are Mott insulators, many of which are driven super- 
conducting by the application of hydrostatic pressure 
or uniaxial stress. Most salts display antiferromagnetic 
(AFM) order, but recent experiments [1, 2] suggest that 
Me3EtP[Pd(dmit)2]2 (P-1) is a valence bond solid (VBS) 
and Me3EtSb[Pd(dmit)2]2 (Sb-1) is a type II spin liquid 
(SL) [2]; with a singlet gap, but no triplet gap [4]. 

In this Letter we report DFT calculations of the elec- 
tronic structures of Sb-1 and P-1. We parameterize these 
results in terms of tight binding models and report the 
parameters found for a number of Pd(dmit)2 salts with 
AFM or charge ordered (CO) ground states. The sim- 
plest model that has been proposed for the insulating 
phases of the Pd(dmit)2 salts is the Heisenberg model on 
the anisotropic triangular lattice [1, 2, 5]. In this model 
each site represents a Pd(dmit)2 dimer, J is the exchange 
coupling around the sides of square and J' is the exchange 
interaction along one diagonal. We find that those mate- 
rials that display long range AFM order lie in the param- 
eter regimes jyj < 0.5 or jyj > 1 where many-body 
theories predict long range magnetic order. Further, all 
of the materials with CO, SL or VBS ground states lie 
in the parameter regime 0.5 < J' / J < 0.9 where the low 
energy physics remains controversial because there are a 
number of competing states. We argue that this means 
that other terms in the Hamiltonian may be crucial for 
determining the ground state. 

The Heisenberg model on the anisotropic triangular 
lattice has been studied by a range of theoretical methods 
including linear spin- wave theory [6], series expansions 
[7], the coupled cluster method [8], large-N expansions 
[9], variational Monte Carlo [10], resonating valence bond 
theory [11-15], pseudo-fermion functional renormaliza- 
tion group [16], slave rotor theory [17], renormalisation 
group [18], and the density matrix renormalisation group 
[19]. Collectively these studies suggest that Neel (7r,7r) 
order is realised for jyj < 0.5 and incommensurate 




FIG. 1. (Color online.) In-plane crystal structures and inter- 
dimer hopping integrals for Sb-1 (left) and P-1 (right). Rings 
indicate the dimers. The conventional unit cell is marked. 



(g, g) long range AFM order is realised for J'/ J ~ 1 (in 
the special case J' / J = 1 the 120° state with q = 27r/3 is 
realised). However, the ground state for 0.6 < J' / J < 0.9 
remains controversial. 

Many Pd(dmit)2 salts undergo a Mott transition under 
hydrostatic pressure and/or uniaxial strain [1, 2]. There- 
fore, it is possible that the Heisenberg model misses some 
essential physics as it represents only the lowest order 
terms in an expansion in t/U^ where t is the hopping in- 
tegral between neighbouring sites and the Hubbard U is 
the effective Coulomb repulsion between two electrons on 
the same site. The next non-trivial term in this expansion 
introduces ring exchange into the Hamiltonian [2, 20]. It 
has been argued argued [20], in the context of the SL 
(BEDT-TTF)2Cu(CN)3 (/^-CN), that ring exchange and 
other higher order terms can drive the breakdown of the 
120° order found on the isotropic triangular lattice be- 
fore t/U becomes large enough for the Mott transition 
to occur. Thus, it is clear that the ratio t/U is vitally 
important in organic charge transfer salts. This has stim- 
ulated several groups to attempt to calculate the U from 
first principles. However, different methods give some- 
what inconsistent results [21-24] and so the ratio t/U is 
not reliably known at present [2]. Nevertheless different 
methods give reasonably consistent trends in the the vari- 
ation of U across trends. Therefore, we have also inves- 



tigated the trend in the U associated with [Pd(dmit)2]2 
dimers. We find very little variation of U across the se- 
ries, which suggests that ring exchange and other higher- 
order effects are not primarily responsible for different 
physics observed in the different Pd(dmit)2 salts. 

A major impediment to comparing the many body the- 
ories described above to experiment has been the lack 
of understanding how changing the cation, i.e., choos- 
ing Pn and n, changes the parameters in the effective 
Hamiltonians of the monomer and dimer models. The 
only previous parameterizations of the band structures of 
the Pd(dmit)2 salts have come from the extended Hiickel 
model (a parameterized tight-binding model). However, 
it has been discovered, e.g., from studies of BEDT-TTF 
salts [21-25], that the extended Hiickel model does not 
provide parameters that are accurate enough for discus- 
sions of the subtle effects of quantum frustration, which 
are at play in both the BEDT-TTF and Pd(dmit)2 salts 
[1, 2]. Therefore, we performed DFT band structure cal- 
culations in Quantum Espresso [26] using the PEE func- 
tional and ultrasoft pseudopotentials with a plane-wave 
cut-off of 25 Ry a 250 Ry integration mesh. Crystal struc- 
tures (Fig. 1) are taken to be those observed by x-ray 
crystallography [27, 28] with only the position of the hy- 
drogen atoms in the cations (which are not visible to 
x-rays) relaxed. 

In molecular acceptors, such as Pd(dmit)2, the Hub- 
bard U associated with a dimer may be written as the 
sum of two terms U = U^^^ — 6U^^\ where U^'^^ 
/•:{()) + E{-2) - 2E{-1) is the U of the dimer in vacuo, 
E{q) is the ground state energy of the dimer with charge 
q and SU'^^^ is the correction due to the polarizable envi- 
ronment of the molecular solid [24]. We performed DFT 
calculations E{q), and hence t/^^^, in SIESTA [29] usmg a 
triplet-zeta plus polarization basis set [24]. However, for 
solids formed from narrow molecules such as Pd(dmit)2 
and BEDT-TTF the accurate calculation SU'^p^ remams 
an challenging problem [24, 30]. Therefore, we report 
U^'^\ which allows us to understand the trends U across 
the systems discussed below. 

Although the in-plane crystal and electronic structures 
of Sb-1 and P-1 are very similar (cf. Fig. 1), there are 
some subtle differences that should be noted before we 
discuss our results in detail. Sb-1 takes the, bi-layer, 
phase, found in many Pd(dmit)2 salts. In the layer shown 
molecules stack face-to- face along a-b, but in alternate 
layers molecules stack in the a+b direction, where a and 
b are the crystallographic axes of the conventional unit 
shown in the Fig 1. The Wigner-Sietz unit cell of Sb-1 is 
half of this size and therefore only contains one dimer per 
layer per unit cell. P-1 also forms a bi-layer structure, 
however, in Pd(dmit)2 molecules stack along the a+c 
axes in both layers. Further, in P-1, the conventional 
and Wigner-Seitz unit cells are identical and contain two 
dimers per layer per unit cell. Thus the Wigner-Seitz 
unit cell of P-1 is approximately twice the size of that of 
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FIG. 2. (Color online.) Band structure of Sb-1 calculated 
from DFT (points) and fit to the dimer tight-binding model 
(curves). The corresponding Fermi surface and the locations 
of high symmetry points are shown in Fig. 5 of the supp. info. 
The values of the dimer tight-binding fit (in meV) for the 
HOMO band are: ti = 37, t2 = 49, ts = 45, t± = -2.1, and 
the parameters for the LUMO band are ti = 7.2, t2 = —15, 
ts -0.2, t± -1.5. 
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FIG. 3. (Color online.) Band structure of P-1 calculated 
from DFT (points) and fit to the dimer tight-binding model 
(curves). The corresponding Fermi surface and the locations 
of high symmetry points are shown in Fig. 6 of the supp. info. 
The values of the dimer tight-binding fit (in meV) for the 
HOMO band are: ti 45, t2 52, H 52, 6.7, b 16, 

and the fit to the LUMO band yields t\ — 6.6, ti — —11, 
t3 -8, t± -8.0, ^ = 0. 

Sb-1 and contains twice as many molecules. 

In Fig. 2 we report the calculated band structure of 
Sb-1. The electronic structure is rather similar to those 
of N-0 and P-0 [31], which are the only other DFT band 
structures for salts of Pd(dmit)2 that we are aware of. In 
all three salts two bands cross the Fermi energy. These 
bands are derived predominately from the antibonding 
combination of highest occupied molecular orbitals (HO- 
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TABLE I. Calculated parameters for the anisotropic triangular lattice for a range of Me4-nEtnPn[Pd(dmit)2]2 (Pn-n) and 
A«-(BEDT-TTF)2Cu(CN)3 (a^-CN). Where two values are given this indicates the small differences because the two t hopping 
integrals are not identical. Values of ti, t2 and ts for all compounds are given in the supplementary information; where we 
also report equivalent results for the local density approximation, which give the same picture as the PBE results. The column 
labelled expt. summarizes the experimentally observed low temperature physics. U^'"^ is the calculated effective Coulomb 
repulsion between to holes on the same dimer in vacuo. The variation is small across all of the Pd(dmit)2 salts and even 
between the Pd(dmit)2 salts and the BEDT-TTF salt. 



MOs). It is not immediately obvious that the HOMOs 
should be partially filled as Pd(dmit)2 is electron accep- 
tor, one might expect a partial occupation of the lowest 
unoccupied molecular orbitals (LUMOs). However, the 
strong hybridisation between the two molecules in each 
dimer (ringed in Fig. 1) means that the bonding combi- 
nation of LUMOs is lower in energy than the antibonding 
combination of HOMOs [2, 31]. There are two bands de- 
rived predominately from the antibonding combination 
of HOMOs and two bands derived predominately from 
the bonding combination of LUMOs because of the two 
dimers per unit cell. 

In Fig. 3 we report the calculated band structure of 
P-1. There are twice as many bands as there are in 
Fig. 2 because the Winger-Seitz unit cell of P-1 contains 
four dimers, rather than two. The charge densities corre- 
sponding to the bands that cross the Fermi energy in both 
Sb-1 and P-1 are reported in the Sup. Info. We also show 
fits to the dimer models in Figs. 2 and 3. The in plane 
tight binding integrals are marked in Fig. 1. Interlayer 
hopping is described by the hopping integral, t±. For 
P-1 we also introduce a parameter which describes the 
orbital energy differences in the HOMOs (LUMOs) due 
to the crystallographically distinct local environments of 
the two dimers per unit cell per layer. Li et al. [32] 
have argued for a quarter filled model where each site is 
a single Pd(dmit)2 molecule. Fits to the monomer mod- 
els are reported in the Sup. Info. Both models provide a 
good description of the DFT results for both compounds. 
However, the monomer fit does not appear significantly 
better than the dimer fit. Therefore we conclude that, at 
the level of band structure, the dimer model is sufficient 
to describe these salts. 

In both materials t2 — ts > ti. Therefore, in or- 
der to make connection with theories of the Heisen- 
berg model on the anisotropic triangular lattice we make 
the identification t = ^(^2 + ^3), = ti. Thus we 



find that f/t = 0.79 for Sb-1 and t'/t = 0.87 for P- 
1. jyj = {f/tf to leading order in t/U. This yields 
jyj = 0.62 for Sb-1 and J' /J = 0.75 for P-1. Both of 
these values are in the regime where the ground state of 
the Heisenberg model remain controversial. 

We also calculated the band structures of a range of 
other salts of Pd(dmit)2. The values of jyj determined 
analogously to those for Sb-1 and P-1 are reported in 
Table 1. Full details of the calculations will be reported 
elsewhere. We also report our parameterization of the 
band structure of the /^-(BEDT-TTF)2Cu(CN)3 (/^-CN), 
which has a SL ground state. Note that the value of 
f /t for A^-CN is in excellent agreement with other esti- 
mates from DFT (0.83 [21] and 0.80 [22]). We also report 
the calculated values of C/^^^ for dimers of Pd(dmit)2 in 
each of these systems in Table 1. There is no significant 
variation in U'^'^^ amongst the Pd(dmit)2 salts and even 
the BEDT-TTF salt has a remarkably similar C/^^). The 
bandwidths W of salts are also remarkably consistent, 
cf. the Supp. Info., with the exception of P-2, which has 
a rather wider band; however, this material is displays 
long-range antiferromagnetism. Therefore, we find no ev- 
idence that ring exchange is the primary determinant of 
the magnetic ground state in these materials. This con- 
clusions is supported by experiment. None of the AFM 
states have been observed to give way to spin-liquid, va- 
lence bond solid or charge order when pressure is to them, 
which drives the materials towards the Mott transition 
by decreasing U /t and hence increasing the relative im- 
portance of ring exchange [1, 2]. However, note that this 
comparison of the relative importance of ring-exchange 
across the series does not imply that ring exchange does 
not play any role in determining the ground state in the 
region where there are multiple competing ground states, 
as we will discuss below. 

It is clear from Table I that those materials which 
display long range magnetic order lie in the parame- 
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ter ranges J' / J < 0.5 or J' / J > 1, while those ma- 
terials which are found experimentally to have CO, SL 
and VBS ground states all lie in the parameter range 
0.5 < J' / J < 0.9. This is precisely the parameter regime 
where the many-body ground state remains controver- 
sial. This suggests that, in this parameter regime, there 
are a number of competing ground states, and that other 
interactions (not included in the anisotropic triangular 
lattice Heisenberg model) may be important for deter- 
mining the ground state. Indeed series expansions cal- 
culations [7] find a number of ground states with very 
similar energies in this parameter regime, supporting our 
contention. Similarly the energies of different phases are 
found to be similar in this regime in coupled cluster cal- 
culations [8]. 

Since we placed this manuscript on the arXiv, Hauke 
[33] has used the compared the tight binding parame- 
ters reported in Table I to his self-consistent spin-wave 
calculations for the Heisenberg model on the completely 
anisotropic triangular with three distinct exchange in- 
teractions. Taking J2/J1 = (^2/^1)^ and J3/J1 = 
(^3/^1)^ he finds that our parameters for the completely 
anisotropic lattice still predict ground states consistent 
with those found experimentally. 

We therefore conclude that the Heisenberg model gives 
a clear prediction of when long-range magnetic order 
win be found in the Pd(dmit)2 salts. For 0.5 > J' / J 
and J' / J > 0.9 the Heisenberg model has a very 
stable ground state with long range magnetic order 
(as witnessed by the good agreement between different 
quantum-many body theories). Thus, perturbations do 
not change the nature of the ground state and long range 
order is realised. But, for 0.5 < J' / J < 0.9 there are a 
number of states that are close in energy (consistent with 
the disagreement in the ground state predicted by differ- 
ent methods). Thus other terms in the Hamiltonian, such 
as intradimer dynamics, ring exchange, elastic forces in 
the crystal or the differences in the crystal structures 
will play an important role in determining which phases 
are realised. In this light it is interesting to note that the 
LOMO bands are within ^ U of the HOMO bands, which 
suggests that they could play an important role. To that 
end we report the tight binding fits of these bands in 
Figs. 2 and 3. This is in marked contrast from proposals 
that proximity to a putative quantum critical point is the 
determining factor in these systems [5, 34]. 
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As well as the Figs, contained in this document ten movies are also included in the supplementary information 
for this Letter. These movies allow for three dimensional viewing of the electronic densities of the bands that cross 
the Fermi levels in P-1 and Sb-1. http: //www.youtube . coin/watch?v=p2ziIq3WLS0&f eature=plcp, http://www. 
youtube . coin/watch?v=LOcpqNZFAJg&f eature=plcp, http: //www.youtube . coin/watch?v=6dJpaNWf Hks&f eature= 
plcp, and http: //www.youtube. coin/watch?v=7_F8iNR3SF0&feature=plcp show the bands from P-1, and http: 
//www.youtube . coin/watch?v=BaN7WV3BNos&feature=plcp shows the density summed over these for bands, http: 
//www.youtube . coin/watch?v=GRXouylwTTg&feature=plcp is a concatenation of all five of movies for P-1. http: 
//youtu.be/-oER-waMTOE and http://youtu.be/4d3UAjR3zMo show the bands from P-1, and http://youtu.be/ 
cFkfOlOcSOc shows the density summed over these for bands. http://youtu.be/92g0xr_fFrk is a concatenation of 
all the of these movies for Sb-1. 




FIG. 1. The inter- monomer hopping integrals for Sb-1. 

In Figs. 5 and 6 we plot the Fermi surfaces of both materials. There are two quasi- two dimensional sheets in the 
Fermi surface of Sb-1. These arise essentially from the bonding and antibonding combinations of the Fermi surfaces in 
the two layers. {t±^ which is responsible for the hybridisation, is a perturbation as it is an order of magnitude smaller 
than the in plane hopping integrals.) Because the Pd(dmit)2 molecules stack along different direction in alternate 
layers, for t± = one would have two uncoupled approximately elliptical Fermi surface sheets; one 'pointing' towards 
the V and the other towards X (cf. Fig. 5 and Fig. 7). It is interesting to note that both DFT and tight binding 
models implicitly assume coherent charge transport between the layers. However, a number of metallic organic charge 
transfer salts do not display coherent interlayer transport [1]. In the metallic phases of these materials, induced by 




FIG. 2. The inter-monomer hopping integrals for P-1. 
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TABLE I. Parameters for the bands derived predominately from the HOMOs of the salts of Pd(dmit)2 from fits to the PEE 
band structures. Note that as there is no symmetry equivalence between any of the hopping integrals we have taken whichever 
of ti, t2 or is the most different to be t' - in any case this leads to topologically equivalent aniostropic triangular lattices. 
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1.17-1.57 


As-2 


AFM 


40 


44 


51 


0.87-1.11 


0.75-1.23 



TABLE IL Parameters for the bands derived predominately from the HOMOs of the salts of Pd(dmit)2 from fits to the LDA 
band structures. 



the application of pressure [2], this could lead to rather interesting phenomenology. For example, because of the 
complicated shape of the first Brillouin zone (FBZ) in Sb-1 the three dimensional Fermi surface occupies less than 
half of the cross-sectional area of the FBZ in the plane (shown in Figs. 5 and 6) despite being half-filled throughout 
the entire FBZ. Therefore, one might expect a dramatic change in the cross-sectional area observed, for example, in 
quantum oscillation experiments if interlayer coherence is lost and a true two dimensional Fermi surface, which is half 
filled in the plane is formed. More subtle restructuring of the Fermi surfaces may also result in either salt from the 
loss of hybridisation between the multiple sheets of Fermi surface, as this is mediated by t^. The doubling of the 
unit cell in P-1 (compared to Sb-1) halves the FBZ. This results in a folding of the Fermi surface - similar to the 
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FIG. 3. (Color online.) Band structure of Sb-1 calculated from DFT (points) and fit to the monomer tight-binding model 
(curves). The location of high symmetry points are marked in Fig. 4 of the main text. The values of the monomer tight-binding 
fit to the HOMO are (in meV): tm 426, tiL 470, t2H 92, t2L 0.6, t^H 1.8, t^L 17, t^H 68, t^L 19, tn 98, 
tbL — 30, t±H — 1.0, t±L = 1.5, v — —427, 5 — 225, where the subscript H indicates parameters for the HOMO band and the 
subscript L indicates parameters for the LUMO band. 
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FIG. 4. (Color online.) Band structure of P-1 calculated from DFT (points) and fit to the monomer tight-binding model 
(curves). The location of high symmetry points are marked in Fig. 4 of the main text. The values of the monomer tight- 
binding fit are (in meV): tm 603, tiL 493, t2H 84, tsL 32, t3H 124, tsL 24, t4H -68, t4L 5.7, 67, 
t^L = 15, toH = -87, toL = 22, t7H = 105, t7L = 13, tsij = 104, tgL = 15, t±H = -6.5, t±L = 8.4, ly = -400, ^ = 300, where 
the subscript H indicates parameters for the HOMO band and the subscript L indicates parameters for the LUMO band. 




FIG. 5. Fermi surfaces of Sb-1 calculated from the dimer model. 



differences between the Fermi surfaces of the f3 and K: phases of BEDT-TTF salts [3] . Otherwise the Fermi surfaces 
of the two materials are remarkably similar. 
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FIG. 6. Fermi surfaces of P-1 calculated from the dimer model. 




FIG. 7. Fermi surface of Sb-1 calculated from the tight binding dimer model for the parameters fit to the DFT calculation, but 
with = 0. Note that now there is no hybridiation between the two bands that cross the Fermi energy. 
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FIG. 8. Electron density for electrons in band 453 of P-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. Protons are not shown. A rotating animation of this plot, which allows for better visualisation of the entire orbital is 
also available on youtube as part of the supplementary information, see main text for link. 
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FIG. 9. Electron density for electrons in band 454 of P-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. Protons are not shown. A rotating animation of this plot, which allows for better visualisation of the entire orbital is 
also available on youtube as part of the supplementary information, see main text for link. 
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FIG. 10. Electron density for electrons in band 455 of P-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. Protons are not shown. A rotating animation of this plot, which allows for better visualisation of the entire orbital is 
also available on youtube as part of the supplementary information, see main text for link. 
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FIG. 11. Electron density for electrons in band 456 of P-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. Protons are not shown. A rotating animation of this plot, which allows for better visualisation of the entire orbital is 
also available on youtube as part of the supplementary information, see main text for link. 




FIG. 12. Structure of P-1. For comparison with Figs. 6-9. 
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FIG. 13. Electron density for electrons in band 227 of Sb-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. A rotating animation of this plot, which allows for better visualisation of the entire orbital is also available on youtube 
as part of the supplementary information, see main text for link. 
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FIG. 14. Electron density for electrons in band 228 of Sb-1. Colour represents phase of the wavefunction. Spheres are atomic 
nuclei. A rotating animation of this plot, which allows for better visualisation of the entire orbital is also available on youtube 
as part of the supplementary information, see main text for link. 




FIG. 15. Structure of Sb-1. For comparison with Figs. 11 and 12. 
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